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Abstract 

High-level quantum-chemical calculations are reported at the MP2 and CCSD(T) 
levels of theory for the equilibrium structure and the harmonic and anharmonic 
force fields of diacetylene, H — C = C — C = C — H. The calculations were performed 
employing Dunning's hierarchy of correlation-consistent basis sets cc-pVXZ, cc- 
pCVXZ, and cc-pwCVXZ, as well as the AN02 basis set of Almlof and Taylor. 
An empirical equilibrium structure based on experimental rotational constants for 
thirteen isotopic species of diacetylene and computed zero-point vibrational correc- 
tions is determined (rf""^: rc-u = 1.0615 A, rc=c = 1.2085 A, rc-c = 1.3727 A) 
and in good agreement with the best theoretical structure (CCSD(T)/cc-pCV5Z: 
rc-H = 1.0617 A, rc=c = 1.2083 A, rc-c = 1.3737 A). In addition, the computed 
fundamental vibrational frequencies are compared with the available experimental 
data and found in satisfactory agreement. 

Key words: Diacetylene, Butadiyne, Rotation-vibration interaction, Anharmonic 
force field, Equilibrium structure 



1 Introduction 



Diacetylene (DiA), H— C=C— C=C— H, the simplest polyacetylene, is of im- 
portance in organic chemical synthesis [1] and is also an abundant astronomical 
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species. In space, it was detected through infrared spectroscopy in the atmo- 
sphere of Saturn's moon Titan [2] and in the protoplanetary nebula CRL618 
[3]. In the laboratory, it has been the subject of numerous spectroscopic studies 
in the infrared and also in the microwave regime (see Refs. [4,5] and references 
therein). Spectroscopic data were also used to evaluate its structural param- 
eters using different approaches. An advanced experimental near-equilibrium 
{r'^)corr structure was reported by Tay et al. [4] making use of experimental 
rotation-vibration data of a total of nine different isotopic species. Addition- 
ally, various quantum-chemical calculations were performed with an empiri- 
cally corrected frozen-core CCSD(T)/cc-pVQZ calculation [6] representing the 
highest theoretical level employed to date. 

The equilibrium structure of DiA is of interest, as it represents one of the sim- 
plest systems with conjugated carbon-carbon multiple bonds. The question is 
how conjugation affects the structure and in particular shortens the central 
carbon-carbon single bond. It is well known that low-level quantum-chemical 
calculations are not reliable in this respect, with Hartree-Fock calculations 
typically underestimating the conjugation and density-functional theory over- 
estimating it. High-level calculations, preferably at the coupled-cluster (CC) 
level [7] , are needed for reliable predictions. To determine the accuracy of the 
CC calculations, it is essential to provide an equilibrium structure based on 
experimental data. 

In a recent Fourier transform microwave spectroscopic study, the rotational 
spectra of four new isotopologs of monodeutero diacetylene were characterized 
[5] raising the number of known isotopic species of DiA to a total of thirteen. 
In the present work, we report an empirical (or semi-experimental) equilibrium 
structure for DiA together with corresponding high-level CC results. The key 
to the empirical structure is the calculation of accurate anharmonic force fields 
which enables the determination of individual zero-point vibrational contribu- 
tions 



for each isotopolog. In Eq. (1), the vibrational correction Ai^g"'^ is given as 
the sum of calculated rotation-vibration interaction constants ai with di as 
the degeneracy being 1 for non-degenerate (stretching) vibrations and 2 for 
doubly degenerate (bending) vibrations. With those corrections, it is possible 
to obtain empirical equilibrium values for the rotational constants according 
to 



The equilibrium structure is then obtained through a least-squares fit of the 
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structural parameters to the estimated equilibrium moments of inertia le (see, 
e.g., Ref. [8]). In case of DiA these are the three bond distances rc-H, "'^c^C; 
and rc-c- 



2 Computational Details 



Quantum-chemical calculations were performed using second-order M0ller- 
Plesset (MP2) perturbation theory as well as the coupled- cluster singles and 
doubles (CCSD) approach augmented by a perturbative treatment of triple 
excitations (CCSD(T)) [9]. The calculations were carried out using basis sets 
from Dunning's hierarchy of correlation-consistent bases [10]. To be more spe- 
cific, the cc-pVXZ (X=D— 6) sets [10] have been used in frozen-core (fc) MP2 
and CCSD(T) calculations, while the cc-pCVXZ (X=T-5) and cc-pwCVXZ 
(X=T,Q) sets [11,12] have been chosen for the all-electron calculations. Ad- 
ditional calculations were also performed with atomic natural orbital (ANO) 
basis sets from Almlof and Taylor [13]. The latter basis set, AN02[^, was 
employed within the fc approximation. 

Equilibrium geometries were optimized using analytic gradient techniques [14]. 
Harmonic, cubic, and semidiagonal quartic force fields were then obtained us- 
ing analytic second- derivative techniques [8,15]. The required third and fourth 
derivatives for the anharmonic force fields were determined by additional nu- 
merical differentiation as described in Refs. [8,16]. The vibration- rotation in- 
teraction constants ai were calculated using formulas given in Ref. [17] based 
on lowest-order rovibrational perturbation theory. Recalculation of the cubic 
force field is avoided by transforming the computed force field for the main 
isotopolog to the corresponding normal coordinate representations of the other 
isotopologs. Since this is not possible for the semidiagonal quartic force field, 
fundamental frequencies, again determined using lowest-order rovibrational 
perturbation theory [17], are only reported for the main isotopic species. 

All calculations have been performed with the Mainz- Austin-Budapest version 
of ACESll [18]; some of the expensive (larger) calculations were made possible 
by using our recent parallel implementation of CC energy as well as first and 
second derivative calculations [19]. 



^ The AN02 set consists of a (13s8p6d4f2g/5s4p3d2flfg) contraction for C and a 
(8s6p4d3f/4s3p2dlf) contraction for H. 
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3 Source of Experimental Data and Least-Squares Fit 

Experimental Bq values of thirteen isotopologues studied through rotation- 
vibration [4,20] and rotational spectroscopy [5] are collected in Table 1. The 
least-squares fits for the structural parameters have been performed including 
the data for all thirteen isotopologues. The fits have been performed with 
respect to the moments of inertia using the same weight for all considered 
isotopic species. 



4 Results and Discussion 

Table 1 summarizes in addition to the experimental rotational constants Bq 
also the corresponding calculated vibrational corrections as obtained at the 
(fc)MP2/cc-pVTZ, (fc)CCSD(T)/cc-pVTZ, (fc)CCSD(T)/cc-pVQZ, CCSD(T)/cc- 
pCVQZ, and (fc)CCSD(T)/AN02 levels of theory. The vibrational corrections 
ABq have been also computed with other basis sets (not given here explicitly). 
Analysis of all these calculations allows one to draw the following conclusions 
(numbers given in the following with respect to the parent isotopic species): 
(a) the use of polarized split- valence basis sets (i.e., cc-pVXZ) in all-electron 
calculations leads to too large corrections (e.g. 3.80 MHz for MP2/cc-pVTZ 
in comparison to 2.1 MHz for (fc)MP2/cc-pVTZ); (b) MP2 calculations over- 
estimate the vibrational corrections in comparison with CCSD(T) computa- 
tions (e.g., 2.1 MHz for (fc)MP2/cc-pVTZ and 0.91 MHz for (fc)CCSD(T)/cc- 
pVTZ); (c) the use of quadruple-zeta basis sets is recommended, as the smaller 
cc-pVTZ and cc-pCVTZ basis sets yield vibrational corrections with rather 
large errors (e.g. 0.91 MHz for (fc)CCSD(T)/cc-pVTZ in comparison with 1.92 
MHz for (fc)CCSD(T)/cc-pVQZ and 2.22 MHz for (fc)CCSD(T)/cc-pCVTZ 
in comparison with 1.84 MHz for (fc)CCSD(T)/cc-pCVQZ); (d) the use of the 
AN02 set yields similar results as the corresponding cc-pVQZ calculations. 

From this analysis we also conclude that our computed vibrational corrections 
(at the CCSD(T)/cc-pCVQZ and (fc)CCSD(T)/cc-pVQZ levels) should have 
an accuracy of about 0.2 to 0.3 MHz. This accuracy is sufficient for obtain- 
ing an accurate equilibrium structure of DiA from the empirical B^ rotational 
constants and, in turn, it is more than adequate for a theoretical prediction 
of ground state rotational constants Bq. For the latter, the errors in the theo- 
retical determination of high-level Bg values due to remaining basis-set errors 
and still missing electron-correlation contributions certainly are larger than 
the zero-point vibrational correction A^g"''^ [21]. 

Table 2 gathers the new extensive set of structural parameters for DiA. There, 
we report the different empirical structures of DiA as obtained from using the 
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vibrational corrections given in Table 1 and the pure theoretical structures 
based on high-level CC calculations using a variety of basis sets. In addition, 
structural parameters from the literature are given [4,5,6]. 

Comparing our five different empirical rf^^ structures, we note that values of 
1.062 A, 1.208 A, and 1.373 A are obtained for the C-H distance, the C=C 
triple bond, and C-C single bond, respectively, irrespective of the used vi- 
brational corrections. The differences in the vibrational corrections affects the 
bond lengths only in the fourth decimal. Nevertheless, the changes when going 
from the cc-pVTZ to the cc-pVQZ or cc-pCVQZ corrections arc not entirely 
negligible and we consider the results obtained with the (fc)CCSD(T)/cc- 
pVQZ, CCSD(T)/cc-pCVQZ, and CCSD(T)/AN02 corrections the most reli- 
able. The high quality of these three sets of equilibrium distances is also seen 
through the fact that the residuals in the fits are as small as 6 kHz for the 
maximum and 0.7 kHz for the root mean-square deviations. The statistical 
uncertainties in all of the fits arc small being 0.0004 A (CC single bond), 
0.0003 A (CC triple bond), and 0.0001 A (CH bond) and hence in all cases 
below 0.001 A. 

Comparison of the derived empirical r,. structures with the pure theoretical 
geometries reveals a good agreement with the results obtained at our best 
quantum-chemical level (CCSD(T)/cc-pCV5Z). With 0.001 A and less the 
discrepancies are in the expected range and probably due to the neglect of 
higher excitations in the CC treatment [22] . A closer analysis of the extensive 
set of theoretical data in Table 2 furthermore demonstrates the importance of 
core correlation, as the (fc)CCSD(T) calculations yield too long bond distances 
(by about 0.0015 A for the CH distances and 0.003 A for the CC distances). 
However, core-correlation effects can be accurately treated using an additivity 
assumption, since the distances obtained via 

re ~ re(pV6Z, fc) + re(pwCVQZ, ae) - re(pwCVQZ, fc) (3) 

are rc-n = 1.0618 A, rc=c = 1.2083 A, and rc-c = 1.3737 A and hence within 
10-4 A of the CCSD(T)/cc-pCV5Z results. 

For the sake of completeness. Table 2 also gives parameters from other equi- 
librium or near-equilibrium structure determinations. Reasonable agreement is 
found between the best empirical equilibrium structures rf^^^ ((fc) CCSD(T) / cc- 
pVQZ, CCSD(T)/cc-pCVQZ, and (fc) CCSD(T)/AN02) and the near-equilibrium 
{''"m)corr structure of Tay et al. [4] with the biggest discrepancy being in the 
central C— C bond. Somewhat larger discrepancies are found in comparison 
against the structure from Ref. [5]. However, the latter structure was only 
based on a rough empirical correction for the difference between the and Tg 
structure of HC3N (without determining the vibrational corrections explicitly) 
and hence is not comparable from a rigorous point of view. Good agreement 
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is found for the structure given by Botschwina and Puzzarini [6] where a 
theoretical structure obtained at the CCSD(T)/cc-pVQZ level of theory was 
empirically corrected through comparison against an empirical equilibrium 
structure of HC3N. Finally, a comparison of the r|"*^ structures against the 
To structure in Table 2 reveals the significance of including vibrational correc- 
tions in structural determinations. In the latter structure the CH distances 
are significantly shorter (by about 0.005 A) whereas the CC single bond is 
longer (by about 0.002 A) than in the r^"^^ structures. A similar conclusion 
also holds for the comparison with the structure reported in Ref. [5] for 
HCCCCD. 

Comparison of the rf^^ parameters of diacetylene with those of similar mole- 
cules reveals that the C-H bond lengths found in diacetylene (1.6016 A, this 
work), fluoroacetylene (1.6014 A, Ref. [6]), methyldiacetylene (1.6013 A, Ref. 
[23]), acetylene (1.0618 A, Ref. [24]), HC3N (1.0623 A, Ref. [25]) and methy- 
lacetylene (1.061 A, Ref. [26]) are practically identical. Clearly, conjugation 
effects in DiA manifest themselves in the C=C bond length: whereas in acety- 
lene a value of 1.2029(1) A is found [24], the C=C triple bond is significantly 
longer in diacetylene, namely 1.2084 A. The latter value is in good agree- 
ment with those found in similar conjugated systems. Recent examples are 
CC distances in the ethynyl group of the substituted diacetylenes HC4F and 
CH3C4H (1.2080 A [6] and 1.2085 A [23], respectively) and also the C=C 
distances in branched species such as czs-Hex-3-ene-l,5-diyne and (Z)-pent-2- 
en-4-ynenitrile (1.208(3) A [27] and 1.207(3) A [28]). Given the prototypical 
C=C equihbrium bond length in ethylene of 1.3305 A [29]) the formal C— C 
single bond length in DiA of 1.3727 A found here is very short (and com- 
parable to monofiuorodiacetylene (1.3729 A) and methyldiacetylene (1.3734 
A) indicative of some double-bond character. Furthermore, conjugation i.e. vr 
electron delocalization is more pronounced in DiA (and HC4F/CH3C4H) than 
in HC3N where the C=C bond is found to be shorter (intermediate between 
acetylene and DiA) and the C— C bond longer compared to DiA (Table 1). 

Several infrared studies of the vibrational fundamentals, overtones and hot 
bands of DiA are found in the literature (see Ref. [4] and references therein) . 
In Table 3 we report harmonic and fundamental frequencies of DiA as obtained 
through CCSD(T) calculations of the cubic and semi-diagonal quartic force 
fields using the cc-pVTZ. cc-pVQZ, and AN02 basis sets in the fc approxi- 
mation and the cc-pCVQZ basis correlating all-electrons. The accuracy of the 
reported fundamental frequencies obtained using Dunning's quadruple-zeta 
basis sets is probably only a couple of wavenumbers, as it turned out impossi- 
ble to converge the required CC second-derivative calculations at the displaced 
points with very tight convergence thresholds. In addition, we also report the 
corresponding (fc)MP2/cc-pVTZ results. Comparison of the theoretical an- 
harmonic (fundamental) frequencies against the experimental values collected 
in Ref. [30] reveals that the results obtained with the cc-pVXZ (X—T, Q) sets 
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turn out less reliable, while better agreement with experiment is seen for the 
other two basis sets. For both, the cc-pCVQZ and the AN02 set, the com- 
puted frequencies match the experimental values within about 10 cm~^. This 
finding again documents the suitability of the ANO basis sets for frequency 
calculations, while the Dunning basis sets (at least when using the valence 
sets) show some deficiencies [31,32,33,34,35]. 



5 Conclusions 

The equilibrium structure of diacetylene has been determined based on the 
combination of experimental rotational constants Bo of thirteen isotopic species 
and zero-point vibrational corrections ABq calculated at various quantum- 
chemical levels. The empirical equilibrium structures obtained agree to within 
10~"^A irrespective of the theoretical level employed. Prom the present study, 
the new recommended equilibrium structure of DiA is tq-u = 1.0615 A, 
rc=c = 1.2085 A, and tq-c = 1.3727 A. This structure is in good agree- 
ment with complementary high-level CC calculations performed here. The 
structural parameters at the highest level of theory (CCSD(T)/cc-pCV5Z) 
are tc-h = 1.0617 A, rc=c = 1.2083 A, rc-c = 1.3737 A. 

Evaluation of cubic and semi-diagonal quartic force fields calculated at the 
CCSD(T)/cc-pCVQZ and (fc)CCSD(T)/AN02 levels yielded harmonic and 
anharmonic vibrational frequencies being in good agreement with experiment, 
while corresponding (fc)MP2 and (fc)CCSD(T) calculations with the cc-pVXZ 
{X—T,Q,) sets are less satisfactory. 
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Table 1 

Experimental rotational constants Bq and computed zero-point vibrational corrections ABq (all values in MHz) for the various isotopo- 
logues of diacetylene. 



ABo 



Isotopic 




(fc)MP2/ 


(fc)CCSD(T)/ 


(fc)CCSD(T) / 


CCSD(T)/ 


(fc)CCSD(T)/ 


Species 


Bn 


cc-pVTZ 


cc-pVTZ 


cc-pVQZ 


cc-pCVQZ 


AN02 


HCCCCH 


4389 3019(391'^ 


2.098 


0.912 


1.922 


1.838 


2.083 


DCCCCD 


3809 2433f66)^ 


-0.207 


—1.110 


—0.203 


—0.237 


—0.005 


jjl3(-.13(-.13(-.13(-;jj 

-1 — L V__y A — L 


4098 8959f36)^ 


2.145 


1.064 


1.970 


1.890 


2.108 


JJ13Q13QCCH 


4243.7325(111)^ 


2.130 


0.997 


1.953 


1.871 


2.103 


H^^CCCCH 


4258 5465(105)^ 


2.112 


0.953 


1.912 


1.829 


2.063 


HC^^CCCH 


4371.6291(45)'' 


2.113 


0.955 


1.962 


1.879 


2.121 


H13C13C13CCH 


4224.7392(99)^ 


2.138 


1.032 


1.987 


1.905 


2.134 


jjISqISqqISqjj 


4115.0556(42)^ 


2.133 


1.027 


1.935 


1.855 


2.075 


HCCCCD 


4084.45342(7)-'' 


0.815 


-0.219 


0.745 


0.688 


0.927 


H^^CCCCD 


3964.11797(17)-^ 


0.871 


-0.141 


0.774 


0.718 


0.946 


HC^^CCCD 


4066.49893(16)-^ 


0.830 


-0.178 


0.782 


0.725 


0.962 


HCC^^CCD 


4071.64202(16)-^ 


0.851 


-0.160 


0.800 


0.743 


0.980 


HCCC^^^CD 


3977.69016(15)^ 


0.870 


-0.144 


0.775 


0.718 


0.947 



" Ref. (fc)CCSD(T)/cc-pVQZ 
" Rcf. (ac)CCSD(T)/cc-pCVQZ 



Ref. [20]. 
" Ref. [4]. 
f Ref. [5]. 



Table 2 



Equilibrium structures of diacetylene and related molecules (A). 



Method 




rc=c 




(fc)MP2/cc-pVTZ 


1.0620 


1.2194 


1.3687 


(fc)CCSD(T)/cc-pVTZ 


1.0638 


1.2149 


1.3789 


(fc)CCSD(T)/cc-pVQZ 


1.0633 


1.2119 


1.3769 


(fc)CCSD(T)/AN02 


1.0631 


1.2118 


1.3766 


(fc)CCSD(T)/cc-pVQZ+corr (Ref. [6])" 


1.0615 


1.2087 


1.3720 


(fc)CCSD(T)/cc-pV5Z 


1.0630 


1.2111 


1.3764 


(fc)CCSD(T)/cc-pV6Z 


1.0630 


1.2109 


1.3762 


CCSD(T)/cc-pCVTZ 


1.0633 


1.2121 


1.3770 


CCSD(T)/cc-pwCVTZ 


1.0630 


1.2111 


1.3763 


CCSD(T)/cc-pCVQZ 


1.0621 


1.2091 


1.3742 


CCSD(T)/cc-pwCVQZ 


1.0620 


1.2089 


1.3741 


(fc)CCSD(T) /cc-pwCVQZ 


1.0632 


1.2115 


1.3766 


CCSD(T)/cc-pCV5Z 


1.0617 


1.2083 


1.3737 


ro 


1.0561 


1.2079 


1.3752 


{r^)corr (Ref. [4]) 


1.0613(1) 


1.2096(1) 


1.3708(2) 


re (Ref. [5]) 


1.0609 


1.2104 


1.3709 


((fc) MP2/cc-pVTZ) 


1.0623 


1.2077 


1.3736 


rr^((fc) CCSD(T)/cc-pVTZ) 


1.0620 


1.2083 


1.3732 


rr^((fc) CCSD(T)/cc-pVQZ) 


1.0616 


1.2084 


1.3727 




i.UOiO 


i.ZUoO 


1 ^797 
1.0 1 Z ( 


r|"'^((fc) CCSD(T)/AN02) 


1.0614 


1.2084 


1.3726 


Acetylene (Ref. [24]) 


1.0618 


1.2029 




rT"P HC3N (Ref. [25]) 


1.0623 


1.2059 


1.3761 


rT"P HC4F (Ref. [6])^ 


1.0614 


1.2080 


1.3731 


rT"P H3C-C4-H (Ref. [23]) 


1.6013(3) 


1.2085(6)/1.2091(16) 


1.3734(14) 


" Empirically corrected ah initio structure 


, see text. 







^ rCsC given refers to the ethynyl group. 
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Table 3 

Vibrational fundamentals of diacetylene (in cm"-*^). 

Experiment'' (fc)MP2/cc-pVTZ (fc)CCSD(T)/cc-pVTZ (fc)CCSD(T)/cc-pVQZ CCSD(T)/cc-pCVQZ (fc)CCSD(T)/AN02 







lidji 111. 


dlillcll Hi. 


lldjl Hi. 


dliiicli iii. 


licxi lii. 


O n rl O TTYl 

dliiicli iii. 


rl O T*TY1 
iidi iii. 


O n rl O T*TY1 

diiiidi iii. 


rl o VTn 
iidi iii. 


o Ti rl o ywi 
diiiidi iii. 




3332 


3486 


3361 


3458 


3329 


3457 


3281 


3463 


3333 


3462 


3330 




2189 


2194 


2147 


2233 


2188 


2235 


2184 


2243 


2197 


2237 


2190 




872 


897 


867 


887 


849 


892 


855 


894 


865 


891 


861 


MK) 


3334 


3482 


3369 


3454 


3339 


3458 


3288 


3465 


3338 


3461 


3329 




2022 


2008 


1970 


2051 


2014 


2057 


2016 


2064 


2031 


2056 


2020 




626 


620 


625 


623 


643 


632 


616 


636 


635 


639 


627 




483 


487 


452 


474 


419 


481 


484 


484 


490 


482 


479 




628 


629 


602 


633 


614 


634 


616 


640 


635 


640 


627 




220 


230 


217 


227 


216 


220 


225 


221 


223 


220 


219 



The accuracy of the reported fundamental frequencies obtained using Dunning's quadruple-zeta basis sets is probably only a couple of 
wavenumbers, as it turned out impossible to converge the required CC second-derivative calculations at the displaced points with very 
tight convergence thresholds. ^ Ref. [30] and references therein. 



